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Using a density-functional based algorithm, the full IR and Raman spectra are calculated for the 
neutral TisCi2 cluster assuming geometries of Th , Td , D2d and Gz v symmetry. The pentagonal 
dodecahedron is found to be dynamically unstable. The calculated properties of the relaxed structure 
having C3„ symmetry are found to be in excellent agreement with experimental gas phase infrared 
results, ionization potential and electron affinity measurements. Consequently, the results presented 
may be used as a reference for further experimental characterization using vibrational spectroscopy. 



I. INTRODUCTION 

In 1992, Castleman and co-workersQ! discovered the 
first of a family of magic number clusters, called met- 
allocarbohedrenes, amidst a metal-carbon cluster distri- 
bution produced using a laser plasma reactor source. 
Metallocarbohedrenes, or Met-Cars for short, are clus- 
ters of stoictdflmetry M 8 Ci 2 (M = V, Zr, Hf, Ti, Nb, 
Mo, Fe, Cr)EM or Tig^M^glM = Zr, Hf, Nb, Mo, 
Ta, W, Si, Y, Nb, Mo, Ta, W)t3&0. Bearing an unusual 
2:3 metal-carbon reduced stoichiometric ratio, a pentag- 
onal dodecahedral cage structure having point group 
symmetry was proposed by Castleman to account for lig- 
and titration experiments which suggested that the eight 
metal atoms-had similar coordination^. Shortly there- 
after, DanceOu proposed a tetracapped tetrahedron of 
Td symmetry and later showed that the dodecahedron 
was higher in energy by 15 eV. Other geometries have also 
been proposed including structures belonging to the C3„ 
point groupEJOE 2 ]^ 3 ! Comparison of the drift time of 
mass selected TigC^ with calculated mobilities for vari- 
ous structures reinforces the proposed hollow cage s , tr , uc .- 
tures instead of more closely packed cubic structurestillij. 

Despite the considerable amount of interest in the 
chemical and physical properties of Met-Cars, the equi- 
librium geometry is still not known! 3 ]. If these MgCi2 
clusters are to be realized as new materials, produc- 
tion and isolation in macroscopic quantities is impera- 
tive. Early experiments by Castleman et al. showed that 
Met-Cars exposed to air were stable fcpshort periods of 
time possibly due to protection by soott3. Most recently, 
Selvan and Pradeep also approached the bulk synthesis 
problem, using an arc-dischargeEZl Infrared spectra of 
the soot after exposure to air over the course of 30 min- 
utes showed the disappearance of a band at 665 cm -1 
and the emergence of others at 1457, 1122, and 872 cm -1 . 
The "appearing" bands were assigned to methylene wag- 
ging, C-C bond stretching, and methylene rocking modes, 
respectively. As discussed below, the highest energy " ap- 
pearing" mode (1457 cm -1 ) also is close to the energy 
range associated with Met-Cars. Time-lapse variations 
across IR spectra cannot be used to declare the presence 
of a "new" material when the disappearing or emerging 



peak(s) are not known or expected to be signatures of 
the species under investigation. 

Recently, the infrared spectra for neutral Tig C 12 , as 
well as for TigCn and Tii4Ci3, have been measured in 
the gas phase by infrared jaultiphoton excitation using 
a pulsed linear acceleratorlia. Though the measurement 
represents the vibrational signature for a free Met-Car, 
it serves as the first available experimental reference ap- 
plicable for use in the identification of the presence of 
TigCi2 in. unpurified samples. Recently, Gueorguiev and 
PacheccO have calculated the infrared absorption spec- 
trum for the Tig C 12 . Their calculations have shown 
the Met-Cars with distorted Td geometry to be the low- 
est energy ones while the experimental IR-absorption 
spectrum compares better with the calculated spectrum 
for the D2d structure. This is clearly a finite temper- 
ature effect. However, the width of the experimental 
peak near 1400 cm -1 is much larger than predicted byj 
the above mentioned article. Neither the experimentaltSI 
or theoreticaltJ works give strong evidence that there 
is Met-Car infrared intensity near the 665 cm" 1 range. 
Therefore an additional calculation aimed at this ques- 
tion is useful. Another motivation of the present work is 
to identify the symmetry of the lowest energy structure 
within the density functional theory. A comparison of the 
theoretical IR absorption spectra with experiment will 
certainly be helpful in this regard. However, due to the 
fact that the existing experimental measurements were 
carried out at finite temperature, the characterization of 
the ground state structure through IR measurement may 
be unreliable. We show in the following that at finite 
temperature in which the experiments are carried out, 
the spectra is a mixture of low lying isomers leading to a 
large width of the high frequency peak. In this paper, we 
present the complete vibrational modes for a T)2d struc- 
ture and a C3„ trigonal pyramid Tig C 12 as determined 
from all-electron density functional studies. Apart from 
the IR absorption spectrum, we also present our calcu- 
lated Raman scattering spectrum for these structures. 
This is the first calculated Raman spectrum of TigCi2 . 
The reliability of our predicted Raman frequencies can 
be assessed based on agreement between the calculated 
IR spectrum and the gas phase IR results or other exper- 
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imentally determined properties. 



II. COMPUTATIONAL DETAILS 

The geometrical and electronic structures of 
Ti§Ci2 molecule of different stmctures were stud- 
ied using density functional theoryEH. Tha,calculations 
were performed with the NRLMOL codeEl within the 
generalized gradient approximation to the exchange- 
correlation energj£3. The Kohn-Sham orbitals were 
expanded in a basis of linear combination of atomic 
orbitals where each orbital is expanded in a Gaussian 
basis set centered on the atoms. An optimized basis 
setE3 was employed to speed up the calculations. The 
basis for the Ti atom consists of 7 s-type, 5 p-type and 
4 d-type orbitals each of which are constructed from 
a linear combination of 19 Gaussians. For C atom, 5 
s-type, 4 p-type and 3 d-type functions were used which 
in turn are composed of 12 Gaussians. A variational 
integration mesh was used for the analytical compu- 
tation of the potentialcj. Spin-polarized calculations 
were performed taking all electrons into account. The 
optimization of the clusters was symme|tejy adapted. The 
massively parallel version of the code=3 was used and 
the calculations were carried out on a cluster of PCs. 
The geometry optimizations were carried out by using 
the limited-memory Broyden-Fletcher-Goldfarb-Shanno 
(LBFGS) scheme of minimization. The geometry opti- 
mizations were carried till the forces on the atoms were 
below 0.001 a.u. For each geometry, the self-consistency 
cycle was carried out till the energy differences were less 
than 0.000001 a.u. 

To calculate the vibrational frequencies, a dynamical 
matrix was constructed by first displacing the atoms by 
± 0.08 a.u. and calculating the forces. The vibrational 
frequencies of the molecules were calculated by a direct 
diagonalization of the dynamical-matrix. The details of 
the method can be found in Ref.cj. The Raman and in- 
frared intensities were obtained from the derivatives of 
the dipole moment and the polarizability tensor which 
in turn were obtained by applying an electric field of 
strength 0.005 a.u.. For a high precision calculation of 
the vibrational frequencies and the Raman and IR in- 
tensities, the energies of the displaced geometries were 
converged to a tolerance of l.OxlO -8 a.u.. 



III. RESULTS AND DISCUSSION 

The earlier reports of the Met-Car structure have 
shown the TigCi2 clusters with a T^ geometry to be lower 
in energy than the ona,with T^ symmetry It may be re- 
called here that Dancea have shown a barrierless reaction 
pathway between T/j and the T^ structures he consid- 
ered. The path was constrainedi-with the D2 symmetry 
On the other hand, Chen et alJi^ have proposed a 
structure to be more strongly bound. To obtain possi- 



ble candidate structures for the low energy isomers, we 
have performed the geometry optimizations for different 
structures with T^ , T/j , Cz v , D 2 d ,D 2 ^ and D 2 point 
group symmetries. Calculations for the lower symmetry 
structures albeit similar to the higher symmetry ones, 
were carried out so as not to rule out any low symme- 
try structures. The lowest energy structures, which are 
unique are shown in Fig. [l[ We have carried out the op- 
timizations for two structures with the T)2d symmetry in 
which the carbon dimers are oriented differently. We use 
a * to distinguish between them. The one labeled as D2d 
has a structure close to the C3,, structure and therefore 
is not shown in Fig. [lj Similarly, we have not shown 
the structure of the D2h geometry since it is close to the 
Th structure. Also, one of the T^ structures, labeled as 
T d is close to the C3U and therefore is not shown. 

The geometry optimization has revealed the structure 
with C^ v symmetry to be the lowest energy structure 
among the structures considered here. The energies of 
all the calculated optimized structures relative to the en- 
ergy of C^ v structure are given in Table |. We wish to 
point out here that a similar structure with T^ symmetry 
(T^) was found to be slightly higher in energy ( 0.13eV). 
The lower order of the C3,, symmetry group removes 
some of the degeneracies at the Fermi level. A fully un- 
symmetrized optimization leads to a structure lower in 
energy by only 0.01 eV. Therefore, we accept the C^ v 
structure as the lowest energy structure among the ones 
examined in this work. We believe that this structure 
is the same as the distorted T^ reported in Ref.E^I .The 
Trf and T/j structures shown in Fig. || are significantly 
higher in energy, respectively, by 14.94 and 14.5 eV than 
the C3,, structure. The energy difference between the Td 
and T^ isomers is small (0.44 eV) . The D2d structure 
is also found to lie higher in energy by 0.06 eV than the 
C3,, structure. The other structures with and D2 
symmetry are still higher by ~ 2 eV. 

To get more information about the stability of the dif- 
ferent conformers, we have calculated the atomization 
energies and the vertical ionization potential (vIP) as 
well as the vertical electron affinities (vEA). These val- 
ues , along with the magnetic moments, and the HOMO- 
LUMO gaps are summarized in Table [|. We found that 
while the T^ structure has magnetic moment 4 fj,s, the 
other structures are found to possess a smaller magnetic 
moment of 2 [Lb- The calculated binding or atomisation 
energy of all the different structures are high which sug- 
gest that these clusters are highly stable. The vertical 
IP and EA are calculated by assuming the structure of 
the charged cluster to be same as the neutral one. The 
relative values of the vIP and vEA indicate the relative 
stability of the clusters with different geometry. Surpris- 
ingly, the higher energy isomers exhibit higher ionization 
potentials. The vertical IP of the T^ , T>2d and the 
C3„ structures range between 4.51-4.61 eV. The earlier 
reported values of calculated vIP for the T^ structure 
ranges between 4.37 - 4.7 eV whereas the adiabatic IP is 
4.43 eVE3. The vIP or IP reported for other structures 
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are considerably higher than the experimental value. For 
the lowest energy structure with C^ v symmetry, we have 
calculated the adiabatic IP. In this case, the geometry of 
the charged cluster was allowed to relax. The calculated 
ionization potential (4.47 eV) is in excellent agreement 
with the recently determined value (4.4 ± 0.02 eV) from 
near threshold photoionization efficiency curves for the 
Ti 8 C 12 El 

The vertical electron affinity of the Ti Met-Cars were 
reported for Td and Th structures. While the experi- 
mental vertical electron affinity is._L.16 ± 0.05 eV, the 
adiabatic affinity of 1.05 ± 0.05 eVHfl is lower. The vEA 
calculated in the present work for the Cs v , T>2d and T^ 
range between 0.89 - 1.08 eV. The vEA for other struc- 
tures are relatively higher while the D 2 /j shows a very low 
vEA (Table I). The calculated adiabatic electron affinity 
of the C3,, structure is 1.00 eV which is in excellent agree- 
ment with the experimental value of 1.05 eVE_l Apart 
from the G^ v anionic structure, we have optimized the 
geometries of the T d and the &2d structures also. The 
anionic T d and Cs v are degenerate with an energy differ- 
ence of 0.004 eV while the D2d anion lies 0.015 ev above. 
The low electron affinity and the high ionization potential 
of the Tig C 12 signifies the low reactivity of this cluster. 

The energy ordering among the isomers and the good 
agreement with experimental IP and EA gives us con- 
fidence that C3t, structure is a likely candidate for the 
lowest energy structure within GGA. 

Another possible way of identifying the structure will 
be through the infrared absorption spectrum. Recently, 
Heijnsbergen et al. Ej have carried out the measure- 
ment of infrared resonance-enhanced multiphoton ioniza- 
tion spectrum (IR-REMPI) of the Ti§Ci2 clusters within 
the frequency range of 400 cm -1 to 1600 cm -1 . The 
IR-REMPI spectrum closely resembles the conventional 
infrared absorption spectrum in peak position and rela- 
tive intensity. This experiment reports a broad peak of 
the IR-REMPI spectrum centered around 1395 cm -1 . A 
comparison of the calculated spectrum with the measured 
spectrum can shed light on the possible candidate struc- 
tures of the highly stable TisCi2 and also in identifying 
the Met-Cars in a mixed environment. The vibrational 
spectrum of the high lying structures with T^, T^, D2 
and T>2h symmetries have several imaginary frequencies 
which further indicates these structures to be highly un- 
stable. We, therefore, concentrate on the analysis of the 
vibrational modes and the IR and Raman spectrum of 
the two lowest energy structures, namely, the C3t, and 
T>2d structures. Pacheco et al£3 have recently reported 
calculated IR absorption spectrum for the Met-Car with 
Td, Th and T)2d structures and unsymmetric structures 
which are Jahn- Teller distorted. They have shown that 
while the Jahn- Teller distorted T^ structure is the low- 
est energy structure, the absorption spectrum of the D 2( z 
structure matches the experimental spectrum best. They 
suggest that at finite temperature, the clusters with T)2d 
symmetry are most abundant. The calculated spectra in 
Ref.__l were broadened by a Gaussian of FWHM of 40 



cm -1 which merges the closely spaced peaks. For the 
sake of completeness, we also calculated -the IR spectra 
for the D2d structure reported in Rcf.l__l. The start- 
ing geometry differs from the other T>2d structure in the 
orientations of the carbon dimcrs. This structure, upon 
optimization, distorts largely from the starting geometry. 
The final structure is shown in Fig. [I] which we refer to as 
DJj d . This structure, although high in energy with respect 
to the C3„, is stable with no imaginary frequency asso- 
ciated with it. However, the IR spectrum shows peaks 
at frequencies lower than in the experimental spectrum - 
one at 1357 cm -1 and a much larger one at 1310 cm -1 . 
Since this structure lies high in the energy scale, we con- 
centrate on the IR spectra of the geometries labeled as 
Cs v and the D2d . 

Although both the C3,; and the T)2d symmetry groups 
are subgroups of the T^ symmetry group, they are inde- 
pendent. The small energy difference of 0.06 eV between 
them suggests the possibility of a symmetry breaking re- 
action path connecting the C3,, and D2d structures. In- 
deed, there exists such a path shown in Fig. |[ The 
plot shows the small energy difference between the two 
structures. However, an estimate of the vibrational fre- 
quencies for such an anharmonic potential ruled out the 
possibility that at low or room temperature the Met-Cars 
can vibrate between the two structures. This fact again 
establishes that at zero temperature, it is likely to con- 
form to the C3,, structure. 

The experimental spectrum was measured between 400 
- 1600 cm -1 range. Therefore, we present our calculated 
spectrum in this frequency range in Fig. 0. The spectra 
shown in Fig. ^ in the upper two panels were broadened 
with a Gaussian of full width half maximum (FWHM) of 
6 cm -1 as well as 40 cm -1 . The high resolution spectrum 
shows two clear peaks at 1393 cm -1 and 1442 cm" 1 for 
the C3„ structure (Fig. |3| lower panel) while the D2d 
structure shows peaks at 1354 cm -1 and 1393 cm -1 . At 
temperature T=0 K, both the structures have strong ab- 
sorption peaks near 1400 cm -1 . The smearing of the 
peaks with a gaussian of FWHM of 40 cm" 1 merges both 
the peaks and brings the spectra closer to the experimen- 
tal one. In the experimental spectra, we detect a shoulder 
near 1364 cm -1 which is clearly reproduced for the C3„ 
spectra but at a higher frequency. Although the shape of 
the IR absorption spectrum is correctly reproduced for 
the calculated C3 U spectra, the spectra is shifted slightly 
towards high energy region. The peaks of the D2d spec- 
tra at 1354 and 1393 cm" 1 have comparable intensities 
and a smearing with Gaussian of FWMH of 40 cm" 1 
does not show a shoulder clearly and the broad peak is 
to the left of the experimentally observed peak. In the 
low frequency region around 500 cm" 1 , the experimental 
spectrum shows peaks at 455 and 520 cm -1 , which are 
correctly reproduced in case of C3,, structure and slightly 
shifted to higher frequency in case D 2( z structure. Based 
on this observation and as well as the fact that the C3„ 
energy is lower than the T>2d , we believe that at tem- 
perature T=0K the IR spectrum will be dominated by 
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that of the C3,, structure. We have estimated the shift in 
the frequency of the highest peak of the C 3 „ structure 
due to the anharmonicity of the potential to be about 10 
cnr 1 . These estimates also showed that the lowest eight 
excitations would encompass a broadening of nearly 30 
cm . 

Given that the experimental IR spectra was associated 
with clusters hot enough to undergo thermionic emission, 
the effect of thermal fluctuations will be strong in the ex- 
perimentally measured IR spectra. Therefore a statisti- 
cally weighted spectra of all the low-lying structure will 
better reproduce the experimental IR absorption spec- 
trum. In Fig. ||, we show the weighted spectra of the C3,, 
and the T>2d structures which are within a range of 0.06 
eV. In the Fig. |J, we show the calculated weighted IR 
absorption for temperatures T=300K, 750K and HOOK. 
At room temperature, the spectra is dominated by the 
structure and weighted spectra has a peak shifted 
towards the high energy region. However, as the temper- 
ature is increased, the contribution from the T)2d struc- 
ture increases and the spectra shows two close peaks for 
750 K. At still higher temperature of HOOK, the peak due 
to the C3,, is reduced. The noticeable feature is that the 
width of the experimental peak is better reproduced in 
the weighted spectra. The contribution from other iso- 
mers will become important at still higher temperature. 
Another noteworthy point is that the low energy region of 
the spectrum is better reproduced at room temperature. 

For the purpose of comparison with any probable fu- 
ture experimental investigation, we also present our cal- 
culated Raman scattering spectra for both C3„ and D 2c ; 
geometries. This is the first calculated Raman spectra 
for TigCi2 . These plots are shown in Fig. ||. The Ra- 
man active modes are seen at lower frequencies which 
arises due to the Ti-Ti stretch mode and also twisting of 
the C-C bonds. These modes will change the volume of 
the cluster and hence the polarizability which will lead 
to Raman activities. The scattering intensities at high 
frequencies are low for both the structures. The spectra 
shown in Fig. H has been convoluted with a Gaussian of 
FWHM 6 cm . A larger value of FWHM will smear out 
most of the fine structures seen in the plots. A few weak 
peaks are seen in the region around 1400 cm -1 where the 
IR activity is most strong. In the Raman spectra of both 



the C3U and D2<j structures, the most prominent peaks 
occur in the region between 100 to 200 cm -1 which are 
nearly similar (Fig. |J). Therefore, Raman spectra in this 
low frequency region will not help in distinguishing the 
ground state structure between C 3lI and T)2d ■ However, 
the spectra between 300 - 600 cm -1 are distinguishable 
as can be seen from Fig. |^. A measurement of the Ra- 
man spectrum can greatly influence the debate over the 
equilibrium geometry of the Met-Cars. 

IV. CONCLUSIONS 

In conclusion, we have carried out extensive density 
functional calculations on the electronic structure and 
the vibrational states of the TigCi2 for different geome- 
tries. We find the geometry with C3,, symmetry to be 
the lowest energy structure. A study of the vibrational 
states show that the calculated IR absorption spectra 
for this geometry compares best with the experimental 
one. We find another structure with T)2d symmetry as 
a competing structure with an energy difference of only 
0.06 eV. We find that a symmetry breaking reaction path 
exists from the T)2d to the G^y structure. However, 
our estimates of the vibrational modes in the potential 
surface between the C3U and T)2d ruled out any possi- 
bility of cold clusters vibrating between the two closely 
placed geometries. We estimate the reaction barrier to 
be about 600K. We point out that the experiments may 
be associated with hot clusters and show that the overall 
width of the high frequency IR spectra supports a high 
temperature mixture of C3,, and T)2d ■ The two peak 
character observed in our calculations could reduce to a 
shoulder structure with small perturbation of IR inten- 
sity and peak positions. We also present our calculated 
Raman spectra for these two low- lying structures. An ex- 
perimental measurement of the Raman spectra may help 
in deciding the ground state geometry of the Ti§Ci2 . 
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FIG. 1: The optimized geometries of the TigCi2 with various 
symmetries. 




FIG. 2: A possible reaction path from the r>2d to the C3,, 
isomer. The x-axis shows the distances along the reaction 
path from the Cg v isomer. 
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FIG. 3: The calculated IR intensities for Tig C 12. The up- 
per panel shows the absorption spectra for D 2 <j geometry and 
the lower one shows the same for geometry. The solid 
lines represent the experimental spectra. The dash-dotted 
line shows the calculated spectra broadened with a gaussian of 
FWHM of 6 cm" 1 while the dashed line shows the one broad- 
ened with 40 cm - . The convoluted curves are not renormal- 
ized. The theoretical intensities are in (D/A) 2 amu _1 . The 
experimental data are scaled down for comparison with the- 
oretical data. 
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FIG. 4: The calculated IR intensities of the Tig C 12 at dif- 
ferent temperatures. The theoretical curves were convoluted 
with Gaussian of FWHM of 40 cm" 1 . 
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FIG. 5: The Raman scattering spectra of the low-lying D 2 d 
and the C3„ structures. The spectra are broadened with a 
Gaussian of 6 cm -1 . 



